Topological spin-current in non-centrosymmetric superconductors 
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We study the spin transport properties of the non-centrosymmetric superconductor with time- 
reversal-symmetry where spin-triplet (p x ±ip !/ )-wave and spin-singlet s-wave pair potential can mix 
each other. We show that when the amplitude of (p x ± ip H )-wave pair potential is larger than that 
of s-wave one, the superconducting state belongs to the topologically nontrivial class analogous to 
the quantum spin Hall system, and the resulting helical edge modes as Andreev bound states are 
topologically protected. We find that the incident angle dependent spin polarized current flows 
through the interface due to the presence of the helical edge modes. With a weak magnetic field, 
also the angle-integrated current is strongly spin polarized. 
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The topological properties of the electronic states have 
been attracting intensive interests in condensed matter 
physics. Especially, it was highlighted by the discovery 
of the quantum Hall system (QHS) showing the accurate 
quantization of the Hall conductance cth which is related 
to the topological integer P, [3] . 

Recently, the concept of the QHS has been generalized 
to the time-reversal (T) symmetric system, i.e., the quan- 
tum spin Hall system (QSHS) 0, 0, S] • QSHS could be 
regarded as the two copies of QHS for up and down spins 
with the opposite chiralities. In generic case, however, 
the mixture of up and down spins occurs due to the spin- 
orbit interaction, which necessitates the new topological 
number to characterize QSHS 0,0). In QSHS, there ex- 
ist the helical edge modes, i.e., the time- reversal pair of 
right- and left-going one-dimensional modes, which has 
been experimentally demonstrated for the quantum well 
of HgTe system by the measurement of the charge con- 
ductance [6|. 

In the field of superconductivity,the chiral p-wave su- 
perconductors such as Sr 2 Ru04 [7] can be considered 
as an analogue of the QHS, and novel phenomena such 
as one-dimensional Maiorana fermions (real fermions) 
modes at the edge |M |9( and the non-Abelian statistics 
of the vortex [3, Il0|| has been proposed there. Beside 
these issues, the non-centrosymmetric (NCS) supercon- 
ductors such as CePtsSi are a central topic [ll|,[l2|]. Also 
the two-dimensional NCS superconductors are expected 
at the interfaces and/or surfaces due to the strong po- 
tential gradient. An interesting example is the supercon- 
ductivity at LaA10 3 /SrTi0 3 interface [l]|. In NCS su- 
perconductors, the spin-orbit interaction comes into play. 
Especially, Frigeri et al. [l2| have shown that (p x ± ip y )- 
pairing state has the highest T c within the triplet-channel 
in CePtsSi. However, the singlet (s-wave) and triplet (p- 
wave) pairings are mixed, and several novel associated 
properties such as the large upper critical field beyond 
the Pauli limit have been focused on 12j. On the other 



hand, the pure (p x ± ip a )-pairing state has been stud- 
ied from the viewpoint of the superconducting analogue 
of QSHS fl3 |. Therefore, it is an important and urgent 
issue to study the spin transport properties of the NCS 
superconductors from the topological viewpoint. 

In this Letter, we study the spin transport properties 



of the non-centrosymmetric (NCS) superconductor [11 1 
with T-symmetry, where (p x ±ip y )-w&ve and spin-singlet 
s-wave pair potential can mix each other. We show that 
when the amplitude of (p x ± ip y )-w&ve pair potential is 
larger than that of s-wave one, the superconducting state 
belongs to the topologically nontrivial class analogous to 
the quantum spin Hall system, and the resulting helical 
edge modes as Andreev bound states(ABS) are topologi- 
cally protected. We study Andreev reflection at low 
energy, which is determined mostly by the helical edge 
modes, and find the incident angle dependent spin po- 
larized current flowing through the interface. When the 
magnetic field is applied, even the angle-integrated cur- 
rent is spin polarized. 

We start with the Hamiltonian of NCS superconductor 



H = 



HQl) 
-A*(-k) 



A(k) 
-H* (-k) 



with H(k) = a + V(k) ■ <x, V(k) = X(xk y - yk x ), 
ik = h 2 k 2 /(2m) - fi. Here, [i, m, & and A denote chemi- 
cal potential, effective mass, Pauli matrices and coupling 
constant of Rashba spin-orbit interaction, respectively 
[13]. The pair potential A(fc) is given by 



A(fc) = [d(k) ■ &]ia y + iip(k)a y 



(1) 



We choose (p x ± ip y )-wave pair for spin-triplet compo- 
nent with d(k) = A p (xk y — yk x )/ \ k \ [13] and s-wave 
one with ip(k) — A s with A p > and A s > 0. The su- 
perconducting gaps Ai = A p + A s and A 2 =| A p — A s | 
open for the two spin-splitted band, respectively, in the 
homogeneous state 161 ] . 
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However, as seen below, surface states are crucially 
influenced by the relative magnitude between A p and A s . 
Let us consider wave function including ABS localized 
at the surface. Consider a two-dimensional semi-infmte 
superconductor on x > where the surface is located at 
x = 0. The corresponding wave function is given by [ijj 

^>s{x) = exp(ik y y)[ciipi exp(iq£ x x) + c 2 ip2 exp(-iq^ x x) 

+c 3 ip 3 exp(iq£ x x) + c 4 ip 4 exp(-iq 2x x)}, (2) 
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Here, fci and &2 are Fermi momentum of the 
small and large magnitude of Fermi surface 

given by —mX/h 2 + J (mX/h 2 ) 2 + 2m[i/h 2 and 

mX/h 2 + \J (mX/h 2 ) 2 + 2m^i/h 2 , respectively. fci(2)a 
denotes the x component of the Fermi momentum 

&1(2), 



with ku2 
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K(2) 
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The wave func- 



tions are given by T ipi = (ui,—ia 1 1 ui,ia 1 1 vi,vi), 
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(u 2 ,ia 2 1 u 2 ,i^a 2 1 v 2 , 
(v 2 ,ia2 1 v 2 ,ija2 1 u 2 ,-jU2), with 7 
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Here we have introduced a.\ — (k\ x ~ ik y )/ki, 
"2 = (&2s - ik y )/k2, &x = {—ki x - ik y )/ki, and 
&2 = {—k 2x — ik y )/k2- E is the quasiparticle energy 
measured from the Fermi energy. 

By postulating ^g{x) — at x — 0, we can determine 
the ABS. The bound state condition can be expressed by 



(A? - E 2 )(A 2 E 2 ) = ]-ji(E 2 + - A, A, 




' i + C 

'[|(01+<fe)] I i \ < f ± 



4>c <| 4>2 |< vr/2, 



(3) 



(4) 



with f < 1, cos 0i = ki x /k\ and cos 02 = k2 X /k 2 . The 
critical angle 4>c is defined as sin -1 (k\Jk 2 ) . For A = 0, 
eq. ([3]) reproduces the previous result [lq . As seen from 
eq. ©, the ABS including zero energy state is only pos- 
sible for I <f>2 I < (f>c an d 7=1) £<e., A p > A s . The 
present ABS is just the edge state, where the localized 
quasiparticle can move along the edge. The energy level 
of the edge state depends crucially on the direction of the 
motion of the quasiparticle. The inner gap edge modes 
are absent for large magnitude of k y , i.e. 4>2- The param- 
eter regime where the edge modes survive is reduced with 
the increase of the magnitude of A. However, as far as 
we concentrate on the perpendicular injection, the edge 



modes survive as the mid gap ABS 18|, [l9( irrespective 
of the strength of A. If we focus on the low energy limit, 
ABS can be written as 

k 2l 



P LA (1 A sA 



My , 



(5) 



with A s < A p for any A with small magnitude of k y . For 
A s > A p , the presnet ABS vanishes since the value of 
right side of eq. (3) becomes negative due to the nega- 
tive sign of 7 for | E |< Ai and \ E \< A 2 . It should be 
remarked that the present ABS do not break the time re- 
versal symmetry, since the edge current carried by each 
Kramers doublet flows in the opposite direction. Thus 
they can be regarded as helical edge modes, where two 
modes are connected to each other by time reversal op- 
eration. 

Now we give an argument why the superconducting 
state with A p > A s has the ABS from the viewpoint of 
Z2 (topological) class We commence with the pure 
{p x ± zpj,)-wave state without the spin-orbit interaction 
A. Spin Chern number Q for the Bogoliubov-de Gennes 
(BdG) Hamiltonian is 2. Turning on A adiabatically, 
which leaves the T-symmetry intact and keeps the gap 
open, one can arrive at the BdG Hamiltonian of interest. 
Upon this adiabatic change of A, the number of the heli- 
cal edge mode pairs does not change. Then we increase 
the magnitude of A s from zero. As far as A p > A s 
is satisfied, the number of helical edge modes does not 
change, since it is a topological number. However, if A s 
exceeds A p , the helical mode disappears. In this regime, 
the topological nature of superconducting state belongs 
to pure s-wave state without A. It is remarkable, just at 
A s = Ap, one of the energy gap of the quasiparticle in the 
bulk closes, where a quantum phase transition occurs. 

Now we turn to the spin transport property governed 
by the ABS in the NCS superconductors [20(. First, we 
point out that the spin Hall effect, i.e., the appearance 
of the spin Hall voltage perpendicular to the supercon- 
ducting current is suppressed by the compressive nature 
of the superconducting state by the factor of (fcp-A m )~ 2 
{kp: Fermi momentum, A m : penetration depth) 21|. In- 
stead, we will show below that the spin transport through 
the junction between the ballistic normal metal at x < 
and NCS superconductor, i.e., (N/SC) junction, can be 
enhanced by the Dopplcr effect at the Andreev reflection. 
We assume an insulating barrier at x — expressed by 
a delta-function potential US(x). The wave function for 
spin a in the normal metal \& jsiix) is given by 



ty N (x) = exp{ik Fy y)[{ip ia 



E 



a<7,pipa P ) exp(ik Fx x) 



E 



bo-,pi>bp exp(-ik Fx x)] 



(6) 



with 7 

T 



T Vfc T =(1,0,0,0), T ^u= T iP H =(0,1,0,0), 
ij} a} =(0,0,1,0), and T Vvi =(0,0,0,1). The corre- 
sponding g(x) is given by eq. (2). The coefficients 
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a a ^ p and b aiP are determined by postulating the boundary 
condition between ^n(x) and ^s(%) given by \I/Ar(0) = 
* s (0), and £^s(0) ~ £^n(0) = 2£g£tf s (Q). We as- 
sume that the Fermi momentum A:^ of normal metal and 
NCS superconductor before putting A are the same. The 
quantities of interest are the angle resolved spin conduc- 
tance fs{4>) and charge conductance fc(4>) defined by 
M 



fsW = [(I « T ,T I 2 " I ou I 2 " I &T.T I 2 + I &u I 2 ) 

2 ,, 2 2 , C0S "ft 

, cos (A 



(I a u | 2 - | a;, T | 2 - | 6 W | 2 + | 6 ijT |< 

/c(0) = [2 + £(|a CT , p | 2 -| V P | 2 )]^. (7) 



Summarizing these features, we can conclude that the 
presence of the helical edge modes in NCS superconduc- 
tor is the origin of the large angle resolved spin current 
through N/NCS superconductor junctions. However, the 
magnitude of the angle averaged normalized spin conduc- 
tance becomes zero since fs(<f>) = — fs(— 4>) ls satisfied. 



p-wave 



where cf> denotes the injection angle measured from the 
normal to the interface. First we consider pure (p x ±ip y )- 
wave state. In Fig. 1, the angle resolved spin conduc- 
tance is plotted as a function of injection angle </> and 
bias voltage V with E = eV. Note here that the k y 
is related to (j> as k y = kpsiiKj). It is remarkable that 
spin conductance has a non zero value although the NCS 
superconductor does not break time reversal symmetry. 
fs{4>) has a peak when the angle <p or k y gives the energy 
E in the energy dispersion of ABS. With this condition, 
the spin-dependent Andreev reflection occurs to result in 
the spin current. Besides this property, we can show that 
fs{<l>) = — fs(— 4>) is satisfied. By changing the sign of 
eV, fs{<P) changes sign as seen in Fig. 1(a). Next, we 
look at the case where s-wave component coexists. We 
can calculate spin current similar to the pure (p x ± ip y )- 
wave case. For A s < A p , where helical edge modes exist, 
fs((j>) has a sharp peak and fs(<t>) = — fs(— </>) is satisfied 
[see Fig. 1(b)]. These features are similar to those of pure 
(p x ± ipj,)-wave case. On the other hand, for A s > A p , 
where the helical edge modes are absent, sharp peaks of 
fs{4>) as shown in Fig. 1 are absent. 

We have checked that there is negligible quantitative 
change by taking A = limit compared to Fig. 1, e.g., 
less than 0.5% change of the peak height. In this limit, 
for pure (p x ± ip y )-wa,ve state, fs{<j>) is given simply as 
follows 

— 8(7^(1 — (Tn) sin 2(f) sin 2ip cos <j> 
| 4(sin 2 (f> — sin 2 <p) + aN[exp(~2itp)(aN — 2) + 2 cos 20] | 2 ' 

for \ E \< A p and f s (<t>) = for | E |> A p with 
sintys = E/Ap Transparency of the interface er/v is given 
by 4 cos 2 0/(4 cos 2 <fi + Z 2 ) with a dimensionless constant 
Z = 2mV /h 2 kp. The magnitude of fs{4>) is largely en- 
hanced at E = ± A p sin <A corresponding to the energy 
dispersion of ABS. The origin of nonzero fs(<P) even with- 
out A is due to the spin-dependent ABS. We have checked 
that even if we take into account the spatial dependence 
of the (p x ±ipj,)-wave pair potential explicitly, the result- 
ing /s(0) does not qualitatively change [201 ] . 
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FIG. 1: (Color online) Angle resolved spin conductance for 
Z = 5. a: eV = 0.1A P , b: eV = -0.1A P and c: eV = 0.6A P 
with A&f = 0.1/j.. (a)pure (p x ± ip y )-wa.ve case with A s = 0. 
(b)A s = 0.3A P . 



Magnetic field offers an opportunity to observe the 
spin current in a much more accessible way, where T- 
symmetry is broken by the shielding current at the in- 
terface. Here we consider the angle averaged normalized 
spin conductance 175 and charge conductance on as a 
function of magnetic field which are given by 2, 22 1 



<7 C 



(8) 



where fNc{4>) denotes the angle resolved charge conduc- 
tance in the normal state with A p = A s = 0. Now 
we consider the magnetic field H applied perpendicu- 
lar to the two-dimensional plane, which induces a shield- 
ing current along the N/NCS superconductor interface. 
When the penetration depth of the NCS superconductor 
is much longer than coherence length, the vector poten- 
tial can be approximated as A(r) = (0, A y (x), 0) with 
A y {x) = —X m Hexp(—x/X m ) with the penetration depth 
A m . Here we consider the situation where the quantiza- 
tion of the Landau level can be neglected. Then qua- 
siclassical approximation becomes available. The ap- 
plied magnetic field shifts the quasiparticle energy E 
in wave function of ^s(x) to E — HA p sin (f>/Ho with 
H = h/{2eir 2 £\ m ) and £ = h 2 k F /{TrmA p ) [23]. For 
typical values of £ ~ lOnm, A m ~ lOOnm, the magnitude 
of Hq is of the order of 0.2Tesla. Here the order of the 
energy of Doppler shift is given by HA p /Hq. Since the 
Zeeman energy is given by [ibH, the order of the energy 
of Doppler shift is kpXm times larger than that of Zee- 
man energy. Thus, we can neglect the Zeeman effect in 
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the present analysis. This is in sharp contrast to QSHS 
where the Zeeman energy is the main effect of H, which 
opens the gap in the helical edge modes and modulates 
the transport properties Q. The enhanced spin current 
due to the Doppler shift is specific to superconducting 
state not realized in QSHS. 
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FIG. 2: (Color online) Angle averaged spin conductance and 
charge conductance as a function of eV with bias voltage V 
with \k F = 0.1/i. a: H = 0, b: H = 0.2H , c: H = -0.2H Q , 
and d: H — OAHo- Curves b and c of the right panel are 
identical. 

As shown in Fig. 1, to discuss topological nature of the 
helical edge modes, it is sufficient to consider pure (p x ± 
ip y )-wa,ve state. In the following, we choose (p x ± ip y )- 
wave case. In Fig. 2, the spin conductance as and charge 
conductance ac normalized by the charge conductance in 
the normal state are plotted. It should be noted that as 
becomes nonzero in the presence of the magnetic field 
H (see curves b, c and d), since fs(<f>) is n0 more odd 
function of due to the imbalance of the helical edge 
modes. For A = limit, the corresponding helical edge 
modes are given by E — A p (l + H/Hq) sin <j> and E = 
— A p (l — H/Hq) sin^. As seen from the curves b and c, 
the sign of as is reversed by changing the direction of the 
applied magnetic field. On the other hand, the resulting 
charge conductance has different features. For H = 0, 
the resulting line shape of ac is the same as that of chiral 
p-wave superconductor (see curve a of right panel) 0, 



171 l24j . As seen from curves b and c of right panel, ac 



does not change with the change of the direction of the 
magnetic field H. Finally, we show in Fig. 3 the zero- 
voltage as and ac ■ as is nearly linearly proportional to 
H. Note that with a small magnetic field H = OAHo ~ 
lOOOe, as is already of the order of 1. Meanwhile, ac is 
almost independent of H. 

In conclusion, we have studied the spin transport 
property of non-centrosymmetric (NCS) superconduc- 
tor from the viewpoint of topology and Andreev bound 
state (ABS). We have found the incident angle depen- 
dent spin polarized current flowing through the inter- 
face. When the weak magnetic field is applied, even 
the angle-integrated current is largely spin polarized. As 
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FIG. 3: (Color online) Angle averaged spin and charge con- 
ductance for eV — as a function of H with XftF = O.lfi. a: 
spin conductance and b: charge conductance. 



the analogy to quantum spin Hall system (QSHS), the 
ABS in NCS superconductor corresponds to the helical 
edge modes in QSHS. The Andreev reflection via helical 
edge modes produces the enhanced spin current specific 
to NCS superconductor. 
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